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ABSTRACT 

We report the discovery of a planet transiting a moderately bright (V = 12.00) G star, with an 
orbital period of 2.788491 ± 0.000025 days. From the transit light curve we determine that the radius 
of the planet is R p = 1.257 ± 0.053 i? Jup . HAT-P-5b has a mass of M p = 1.06 ± 0.11 Af Jup , similar to 
the average mass of previously-known transiting exoplanets, and a density of p p — 0.66± 0.11 gcm~ 3 . 
We find that the center of transit is T c = 2,454,241.77663 ± 0.00022 (HJD), and the total transit 
duration is 0.1217 ± 0.0012 days. 

Subject headings: stars: individual: GSC 02634-01087 - planetary systems: individual: HAT-P-5b 



1. INTRODUCTION 

To date about 20 extrasolar planets have been found 
which transit their parent stars and thus yield values 
for their mass and radius 10 . Masses range from 0.07 Mj 
(GJ436: iGillon et al.l l2007h to about 9Mj (HAT-P-2b; 
iBakoset al.l [2007l), an d radii from 0.7 flj (G J436) to 
about 1.7.Rj (TRES-4: iMandushev et~aT1l2007t ). These 
data provide an opportunity to compare observations 
with theoretical models of planetary structure across a 
wide range of parameters, in c luding those of the host star 
(e.g. iBurrows et al. I l2007t iFortnev. Marlev. fc Barnesl 
120071 and references therein). Transits also yield pre- 
cise determination of other physical parameters of the 
extrasolar planets, for instance the surface gravity. Inter- 
esting correlations between these parameters have been 
note d early on, such as that betw een masses and peri- 
ods (|Mazeh. Zu ckcr. & Pont 200 5]) or periods an d surface 
gravities (|Southwor th. Whc atlev. fe Samsl[2007l ). Classes 
of these close-in planets have also been suggested, such as 
very hot Jupiters (VH Js; P=l-3 days) and hot Jupiters 
(HJs ; P=3-9 days; iGaudi. Seager. fc Mallen-Ornelasl 
120051 ) , or a possible dichotomy based on Safronov num- 
bers (|Hansen fe Barmanll2007f ). However, the small en- 
semble of transiting exoplanets (TEPs) does not allow 
robust conclusions, thus the addition of new discoveries 
is valuable. 

Over the past year the HATNet project 11 (jBakos et al.1 
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Fig. 1. — The top panel shows the unbinned HATNet light curve 
with 4940 data points, phased with the period P = 2.788491 d, and 
with the binned data overplotted. The 0.013 mag deep transit is 
detected with a dip significance of 18. The other panels show pho- 
tometry follow-up in the following order: Sloan z-band photometry 
taken with the FLWO 1.2 m telescope (on two separate dates: UT 
2007 May 18 and UT 2007 May 21), Cousins R-band photometry 
taken with the Wise 1 m telescope (on 4 nights: UT 2007 May 26, 
June 20, July 4 and 18). Over-plotted are our best analytic fits as 
described in the text. 



120021 l2004f ). a wide-angle photometric survey, has an- 
nounced four TEPs. In this Letter we report on the 
detection of a new transiting exoplanet, which we label 
HAT-P-5b, and our determination of its parameters, such 
as mass, radius, density and surface gravity. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Detection of the transit in the HATNet data 

GSC 02634-01087, also known as 2MASS 
J18173731+3637170 is a G star with / w 11.3 and 
V w 12.00. It was initially identified as a transit 
candidate in our internally labeled field G196, centered 
at a = 18 h 08 h , 8 = 37°30'. The data were acquired 
by HATNet's HAT- 7 telescope at the Fred Lawrence 
Whipple Observatory (FLWO) of the Smithsonian 
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Astrophysical Observatory (SAO) and HAT-9 telescope 
at the Submillimeter Array (SMA) site atop Mauna Kea, 
Hawaii. Following the standard calibration procedure of 
th e frames, data were reduced using the astrometry code 
of iPal fe Bakosl (|2006l ). and a highly fine-tuned aperture 
photometry. We applied our external parameter decor- 
relation (EPD) technique on the light curves, whereby 
deviations from the median were cross-correlated with 
a number of "external parameters" , such as the X and 
Y sub-pixel position, FWHM, hour-angle, and zenith 
distance. We have also applied the Trend Fil tering Algo- 
rithm (TFA; IKovacs. Bakos. fc Novell l2005h along with 
the B ox Least Squares (BLS; IKovacs. Zucker. fc Mazehl 
l2002f) transit-search algorithm in our analysis. For field 
G196 we gathered ~3,750 (HAT-7) plus -890 (HAT-9) 
data-points at 5.5 min cadence between 2005 June 8 
and 2005 December 5 (UT). In the light curve of star 
GSC 02634-01087 we detected a ~13mmag transit 
with a 2.7881 d period, signal-to-noise ratio of 12 in 
the BLS frequency spec trum, and dip-significance of 18 
(|Kovacs fc Bakos 1 12005). The top panel of Fig. Q] shows 
the unbinned light curve with all ~4640 data points, 
folded with the period that we derived subsequently, 
based on high precision follow-up, as described later in 

§EH 

2.2. Early spectroscopy follow-up 

Initial follow-up observ ations were m ade with the CfA 
Digital Speedometer (DS: lLathamlll992l ) in order to char- 
acterize the host star and to reject obvious astrophysi- 
cal false-positive scenarios that mimic planetary transits. 
The four radial velocity (RV) measurements obtained 
over an interval of 33 days showed an rms residual of 
0.41kms _1 , consistent with no detectable RV variation. 
Atmospheric parameters for the star (effective temper- 
ature T e ff, surface gravity log <?, metallicity [Fe/H], and 
projected r otational velocity v sint) were derived as de- 
scribed by iTorres. Neuhauser fc Guentherl (|2002D . The 
first three quantities are strongly correlated and difficult 
to determine simultaneously. For example, the uncon- 
strained value log g = 4.0 ±0.2 we obtained is somewhat 
lower than derived from our stellar evolution modeling in 
§0 which is logg = 4.37. Consequently, in a second iter- 
ation we held log g fixed at this value and redetermined 
the other quantities, obtaining T c g = 5960 ± 100 K, 
[Fe/H] = +0.24 ±0.15, and vsini = 2.6 ± 1.5 kms" 1 ! 
These correspond to a slowly-rotating early G main se- 
quence star. 

2.3. High-precision spectroscopy follow-up 

High-resolution spectroscopic follow-up was car- 
ried out at the Haute Provence Observatory (OHP) 
1.93-m telescope, with the SOPHIE spectrograph 
(|Bouchv fc the Sophie Teanll2006l ). SOPHIE is a multi- 
order echelle spectrograph fed through two fibers, one 
of which is used for starlight and the other for sky 
background or a wavelength calibration lamp. The in- 
strument is entirely computer-controlled and a standard 
data reduction pipeline automatically processes the data 
upon CCD readout. RVs are calculated by numerical 
cross-correlation with a high resolution observed spectral 
template of a G2 star. Similar spectroscopic follow-up 
with SOPHIE has already resulted in t he confirmation 
of two TEPs: WASP-lb and WASP-2b (jCameron et all 



TABLE 1 
Radial Velocities for HAT-P-5. 



BJD -2,400,000 


RV a 




BS b 


(days) 


(ms — - 1 ) 


(ms _1 ) 


(ms- 1 ) 


54227.5199 


7721.4 


12.2 


22.5 


54228.5949 


7457.4 


22.3 


18.2 


54229.6098 


7710.4 


17.2 


-7.0 


54230.4900 


7603.8 


22.4 


-20.0 


54231.6088 


7521.3 


14.3 


15.2 


54233.6057 


7579.4 


15.1 


-31.0 


54234.5210 


7510.5 


21.3 


-54.8 


54255.5171 


7680.5 


9.8 


4.2 



a The RVs include the barycentric correction. Bisec- 
tor spans. 

120071 ). HAT-P-5 was observed with SOPHIE in the high- 
efficiency mode (R ~ 39000) during our May 2-13, 2007 
observing run, with an additional measurement taken on 
June 4. Depending on observing conditions, exposure 
times were in the range of 15 to 35 minutes, resulting in 
signal to noise ratios of 20-55 per pixel at A = 5500 A. 
Using the empirical relation of Cameron et al. (2007) 
we estimated the RV photon-noise uncertainties to be 
10-25ms _1 . We present here 8 radial velocity measure- 
ments taken when the planet was out of transit, listed in 
Table [T] Measurements taken durin g transit, revealin g 
the Rossiter-McLaughlin effect (e.g. IWinn et alj T2005) . 
will be presented in a separate paper. 

2.4. Photometry follow-up 

In order to better characterize the transit parameters 
and also to derive a better ephemeris, we performed 
follow-up photometric observations with 1-m class tele- 
scopes. A partial transit of HAT-P-5b was observed us- 
ing the K eplerCam detector on the FLWO 1.2 m tele- 
scope rsee lHolman et alj|2007l ) on UT 2007 May 18. We 
refer to this event as having transit number N tr = — 1. 
Three days later a full transit, N tI — 0, was observed 
with the same instrument. The two Sloan z-band light 
curves are shown in the middle panel of Fig. [TJ We 
also gathered data for four subsequent full transit events, 
N tI = 2, N t r = 11, N tT = 16 and N tI = 21, using 
the Wise 1 m telescope in the Cousins R band (lower 
panel of Fig. []}. Data were reduced in a similar man- 
ner to the HATNet data, using aperture photometry 
and an ensemble of ~300 comparison stars in the field. 
Since the follow-up observations span 22 transit cycles 
(~2 month time-span), we were able to obtain an ac- 
curate ephemeris. An analytic model was fit to these 
data, as described below in § [SJ and yielded a period 
of 2.788491 ± 0.000025 d and reference epoch of mid- 
transit T c = 2,454,241.77663 ± 0.00022 d (HJD). The 
length of the transit as determined from this joint fit 
is 0.1217 ± 0.0012d (2 hours, 55 minutes), the length 
of ingress is 0.0145 ± 0.0007d (20.9 minutes), and the 
central transit depth is 0.0136 mag. 

3. STELLAR PARAMETERS 

The mass (M p ) and radius (Rp) of a transiting planet 
scale with those of the parent star. In order to determine 
the stellar properties needed to place M p and R p on an 
absolute scale, we made use of stellar evolution mod- 
els along with the observational constraints from spec- 
troscopy. Because of its relative faintness, the host star 
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TABLE 2 

Summary of stellar parameters for 



HAT-P-5. 


Parameter 


Value 


Source 


Teff(K) 


5960 ± 100 


DS 


■usini(kms 1 ) 


2.6± 1.5 


DS 


log a 


4.368 ± 0.028 


Yonsci-Yale 


[Fe/H](dex) 


+0.24 ±0.15 


Yonsei-Yale 


Distance (pc) a 


340 ± 30 


Yonsei-Yale 


Mass (M ) 


1.160 ±0.062 


Yonsei-Yale 


Radius (Rq) 


1.167 ±0.049 


Yonsei-Yale 


log(i*/i ) 


0.187 ±0.064 


Yonsei-Yale 


My 


4.32 ±0.18 


Yonsei-Yale 


Age (Gyr) 


2.6 ± 1.8 


Yonsei-Yale 



a Assuming no extinction due to the proximity of 
the star. 

does not have a parallax measurement from Hipparcos , 
and thus a direct estimate of the absolute magnitude is 
not available for use as a constraint. An alternative ap- 
proach is to use the surface gravity of the star, which 
is a sensitive measure of the evolutionary state of the 
star and therefore has a very strong influence on the ra- 
dius. However, logg is a notoriously difficult quantity 
to measure spectroscopically and is often strongly corre- 
lated with other spectros copic paramete r s (see § I2.2[) . It 
has been pointed out by iSozzetti et alj (|2007f) that the 
normalized separation of the planet, a/J?+, can provide 
a much better constraint for stellar parameter determi- 
nation than log g. The a/ quantity can be determined 
directly from the photometric observations with no ad- 
ditional assumptions (other than limb-darkening, which 
is a second-order effect), and it is related to the density 
of the central star. As discussed later in § [BJ an analytic 
fit to the light curve yields a/R± — 7.50 ± 0.19. 

This value, along with T c g and [Fc/H] from §[6l was 
co mpared w i th th e Yonsei-Y ale stellar evolution m odels 
of lYi et all (|200lf) following ISozzetti et all pOOl . As 
described earlier, the initial temperature and metallicity 
from our DS spectroscopy were subsequently improved 
by applying the logg constraint from the models, and 
the isochrone comparison was repeated. This resulted 
in final values for the stellar mass and radius of = 
1.160 ± 0.062M Q and R* = 1.167 ± 0.049 R Q , and an 
estimated age of 2.6 ±1.8 Gyr. We summarize these and 
other properties in Tabled 

4. SPECTROSCOPIC ORBITAL SOLUTION 

Our 8 RV measurements from SOPHIE were fitted 
with a Keplerian orbit model solving for the velocity 
semi-amplitude K and the center-of-mass velocity 7, 
holding the period and transit epoch fixed at the well- 
determined values from photometry. The eccentricity 
was initially set to zero. The resulting rms residual of 
~23.7 m s _1 is somewhat larger than expected from the 
internal errors, and we find that a reduced x 2 value of 
unity necessitates the addition of uncorrelated noise of 
14.4 m s _1 in quadrature, which we attribute to "stel- 
lar jitter". T his level of jitter is consisten t with the 
predictions of iSaar. Butler. &; Marcvl (|1998l ) for a pro- 
jected rotational velocity such as what we measure for 
the parent star. The final fit, with the internal errors 
increased as described above, yields K = 138 ± 14m s _1 
and 7 = 7613.8 ± 9.1ms" 1 . The observations and fitted 
RV curve are displayed in the top panel of Fig. [2l with 
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Fig. 2. — The top panel shows the RV measurements phased with 
the period of P = 2.788491 d. The zero-point in phase corresponds 
to the epoch of mid-transit. Overlaid is the best fit, assuming 
14.4 ms — 1 stellar jitter. The middle panel shows the residuals from 
the fit. The bottom panel displays the line bisector spans on the 
same scale as the upper panel. No variation in the line bisectors is 
seen concomitant with that in the RVs, essentially confirming the 
planetary nature of the transiting object. 

the residuals shown in the middle panel. 

As a test we allowed for the possibility of an eccentric 
orbit and solved for the two additional quantities e cos u 
and esinw, but the results were insignificantly different 
from zero. 

5. EXCLUDING BLEND SCENARIOS 

We have tested the reality of the velocity variations by 
carefully examining the spectral line bisectors of the star 
using our OHP data. If the velocity changes measured 
are due only to distortions in the line profiles arising from 
contamination of the spectrum by the presence of a bi- 
nary with a period of 2.79 days, we would expect the 
bisector spans (which measure line asymmetry) to vary 
with this period and with an amplitude similar to the ve - 
locities (see, e.g.. lQueloz et al.ll200lUTorres et ai] [2005). 
As shown in the lower panel of Fig. [51 the changes in 
the bisector spans are of the same order as the residual 
RV variations, and much smaller than the radial veloc- 
ity semi-amplitude itself. This analysis shows that the 
orbiting body is a planet and rules out a possible blend 
scenario. 

6. PLANETARY PARAMETERS 

The light-curve parameters of HAT-P-5b were deter- 
mined from a joint fit based on the 6 distinct transit 
events, observed with the FLWO 1.2 m and Wise 1 m 
telescopes. A circular orbit was assumed, based on our 
analysis above. We adopted a quadratic limb-darkening 
law for the star, and took the appropriate coefficients 
from IClaretJ (|2004l ). for both the Sloan z and Cousins R 
bands. The drop in flux i n the light curve s was mod- 
eled with the formalism of iMandel fc Ago] (|2002D . us- 
ing the equations for the general case (i.e. not the small 
planet approximation). The adjusted parameters in the 
fit were i) the mid-transit times of the first full transit 
(N tr = 0, T c0 ), and the last full transit (N tr = 21, T c2 i), 
(this is equivalent to fitting for an epoch E and a period 
P), ii) the relative planetary radius, p = R p / R+; iii) the 
square of the normalized impact parameter, b 2 ; iv) the 
quantity Q/R* = tt/.R*(27r/P)/vl — b 2 . From simple ge- 
ometric considerations (,/ R+ and b 2 have an uncorrelated 
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TABLE 3 

Orbital fit and planetary parameters for the 
hat-p-5b system. 



Parameter 


Value 


Period (d) a 


2.788491 ± 0.000025 


T mid (HJD) a 


2,454,241.77663 ± 0.00022 


Transit duration (day) 


0.1217 ±0.0012 


Ingress duration (day) 


0.0145 ±0.0007 


Stellar jitter (ms _1 ) b 


14.4 


7 (ms _1 ) c 


7613.8 ±9.1 


K (ms- 1 ) 


138 ± 14 


e a 







7.50 ± 0.19 


Rpl R* 


0.1106 ± 0.0006 


a Icl (AU) 


0.04075 ± 0.00076 


ip(deg) 


86? 75 ±0? 44 


M p (Mj) 


1.06 ±0.11 


Rp(Rj) 


1.26 ±0.05 


p p (gcm" 3 ) 


0.66 ±0.11 


g p {ms- 2 ) d 


16.5 ± 1.9 


AT C , , N tr = (10~ b d) 


6 ± 27 


AT Ci0 , N tr = 2 (10" 5 d) 


79 ± 58 


AT C , , Ntr = 11 (10" 5 d) 


6 ± 62 


AT C , , N tr = 16 (lCr 5 d) 


-112 ± 84 


AT C , , N tr = 21 (10~ 5 d) 


11 ± 57 



a Fixed in the orbital fit. b Adopted (see text). c The 
7 velocity is on an absolute reference frame. Based 
on only directly observable quantities, sec 
ISouthworth. Wheatlev. fc Sams] J2007D . 

a posteriori probability distribution in parameter space. 
This amounts to an orthogonalization of the fitted pa- 
ra meters, simil a r (alb eit simpler) to the one employed 
bv lBurke et~aLl (|2007l ) for the case of XO-2b. 

We us ed the Markov Cha in Monte Carlo algorithm 
(see, e.g. iHolman et al.l l2007h to derive the best fit pa- 
rameters. Uncertainties were estimated using synthetic 
data sets, by added Gaussian noise to the fitted curve 
at the dates of our observations, and re-solving the light 
curve analytic fit. The magnitude of the noise was taken 
from the white- and red-noise estimations based on the 
real residuals. This process was repeated 1.5 • 10 5 times, 
yielding a good representation of the a posteriori distri- 
bution of the best-fit parameter values. We found this 
method of error estimation to be robust, since it is not 
sensitive to the number of out-of-transit (OOT) points. 

The result for the radius ratio is R p /R± = 0.1106 ± 
0.0006, and the normalized separation is a/R+ = 7.50 ± 
0.19. We found that the a posteriori distribution of b 2 is 
consistent with a symmetric Gaussian distribution, and 
yields b 2 = 0.181 ± 0.040, therefore the orbit is inclined. 
From the inclination, the mass of the star (Table [2]), and 
the orbital parameters (§ 2]) , the other planetary param- 
eters (such as mass, radius) are derived in a straightfor- 
ward way, and are summarized in Table [31 We note that 
a/R+, as derived from the light curve fit, is an important 
constraint in the stellar parameter determination (§13]), 
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which in turn defines the limb-darkening coefficients that 
are used in the light curve analytic fit. Thus, after the 
initial analytic fit to the light curve and the stellar pa- 
rameter determination, we performed another iteration 
in the light curve fit. We found that the change in pa- 
rameters was imperceptible. 

The possibility of transit time variations (TTVs) was 
checked by fitting the center of the transit of the five 
full transit events independently. We found no sign of 
TTV, as the transit times differ by less than 1-cr from 
the expected values (listed in Table [3j . 

7. CONCLUSIONS 

HAT-P-5b is an ordinary hot Jupiter (P = 2.788 days) 
with slightly inflated radius (R p — 1.26i?j) for its mass 
of 1.06Mj, orbiting a slightly metal rich solar-like star. 
The ~20% radius inflation is what current models predict 
for a planet with equilibrium temperature of ^1500K 
(|Burrows et al. Il2007t [Fortnev. Marlev. & Barnesll2007t ). 

HAT-P-5b is more massive than any of the known 
TEPs with similar period (2.5 < P < 3d), such as XO- 
2b, WASP-1, HAT-P-3b, TRES-1, and HAT-P-4b, with 
the exception of TRES-2. The latter is fairly similar in 
mass, radius, orbital period, and stellar effective temper- 
ature. 

However, HAT-P-5b is interesting in that it falls be- 
tween Class I and II, a s defined by the Safronov number 
and T eq of the planet ([Hansen fe Barman! [20071 ). HAT- 
P-5b has a Safronov number of 0.059±0.005 , while Class 
I is defined as 0.070 ± 0.01, especially at T eq ~ 1500K. 
It seems that the additional discovery and characteriza- 
tion of transiting planets of Jupiter and higher masses 
would be very helpful in order to understand these new 
correlations and their reality. 
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